A Finite Difference Domain DecompositionMethod Using Local Corrections for theSolution of Poisson’s Equation
نویسندگان
چکیده
We present a domain decomposition method for computing finite difference solutions to the Poisson equation with infinite domain boundary conditions. Our method is a finite difference analogue of Anderson’s Method of Local Corrections. The solution is computed in three steps. First, fine-grid solutions are computed in parallel using infinite domain boundary conditions on each subdomain. Second, information is transferred globally through a coarse-grid representation of the charge, and a global coarse-grid solution is found. Third, a fine-grid solution is computed on each subdomain using boundary conditions set with the global coarse solution, corrected locally with fine-grid information from nearby subdomains. There are three important features of our algorithm. First, our method requires only a single iteration between the local fine-grid solutions and the global coarse representation. Second, the error introduced by the domain decomposition is small relative to the solution error obtained in a single-grid calculation. Third, the computed solution is second-order accurate and only weakly dependent on the coarse-grid spacing and the number of subdomains. As a result of these features, we are able to compute accurate solutions in parallel with a much smaller ratio of communication to computation than more traditional domain decomposition methods. We present results to verify the overall accuracy, confirm the small communication costs, and demonstrate the parallel scalability of the method. c © 2002 Elsevier Science (USA)
منابع مشابه
The new implicit finite difference method for the solution of time fractional advection-dispersion equation
In this paper, a numerical solution of time fractional advection-dispersion equations are presented.The new implicit nite dierence methods for solving these equations are studied. We examinepractical numerical methods to solve a class of initial-boundary value fractional partial dierentialequations with variable coecients on a nite domain. Stability, consistency, and (therefore) convergenceof t...
متن کاملRESIP2DMODE: A MATLAB-Based 2D Resistivity and Induced Polarization Forward Modeling Software
Forward modeling is an integral part of every geophysical modeling resulting in the numerical simulation of responses for a given physical property model. This Forward procedure is helpful in geophysics both as a means to interpret data in a research setting and as a means to enhance physical understanding in an educational setting. Calculation of resistivity and induced polarization forward re...
متن کاملThomas–fermi and Poisson Modeling of Gate Electrostatics in Graphene Nanoribbon
We describe a simple graphene nanoribbon and bottom gate system and present numerical algorithms for solving Poisson’s and Thomas–Fermi equations for electrons in the graphene nanoribbon. The Poisson’s equation is solved using finite difference and finite element methods. Using the Poisson and Thomas–Fermi equations we calculate an electrostatic potential and surface electron density in the gra...
متن کاملA New Implicit Finite Difference Method for Solving Time Fractional Diffusion Equation
In this paper, a time fractional diffusion equation on a finite domain is con- sidered. The time fractional diffusion equation is obtained from the standard diffusion equation by replacing the first order time derivative by a fractional derivative of order 0 < a< 1 (in the Riemann-Liovill or Caputo sence). In equation that we consider the time fractional derivative is in...
متن کامل